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ABSTRACT 


Loran-C time difference (TD) readings were’ taken 
consecutively at several survey stations near Monterey, CA, 
over a period of three days. One station, Range-7, was 
designated the "known" point of a differential Loran-C 
system. Readings from the known point were used to correct 
readings from a second survey station, Luces Point. A method 
of improving the precision of Loran-C TD readings based on 
the redundancy and relative accuracy of three LOP's was 
developed and applied to the data. Since only one receiver 
waS available, a linear regression of TD vs time _ was 
Calculated and used in the differential correction. 

Based on 496 sets of data taken at 5-second intervals at 
Range-7 in two groups, before and after about 250 readings at 
Luces Point, the absolute accuracy of Luces Point data was 
improved from about 385 m to about 48 m compared to the known 
position of the point. Precision was improved from about 
14.9 m to about 12.6 m circular error (CEP) using the three- 
station correction. Further improvement would probably have 
resulted if two receivers were available for real time 


corrections. 
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I. INTROOOCTION 


A. Purpose 


The purpose of this’ study was to investigate the potential 
accuracy of the Loran-C navigation system using a differential 
correction technique and a correction for the use of three secondary 
stations rather than the usual two. 

Sources of error were first examined and their magnitudes 
estimated. Methods of expressing the accuracy of a determination of 
geographic position were also considered. 

It was determined that position accuracy is composed of two 
factors: the absolute accuracy of the position (which can be 
considered the bias, or systematic error, in the system) amd the 
precision of the individual readings. 

Over longer periods of data collection it was considered possible 
to subtract out the system bias by means of a quasi-differential 
correction based on intermittent time difference readings at a fixed 
point, designated the "known" point of the differential system, and 
to improve the precision of TD readings by taking advantage of the 
redundancy of the reception of three master-secondary pairs. 

For this purpose, a theory of corrections was developed and an 
experiment was designed to determine the effectiveness of the 
corrections. The experimental data was processed by computer, and the 
relative accuracies were calculated. Computer programs for processing 
and displaying data were written for this purpose. 


J, 


B. Method 


To collect data for processing, several known survey points with 
published positions were chosen as locations for collecting Loran-C 
data. Using known points provided a means of comparing results. One 
point in a relatively central location was chosen as the "known" point, 
where a larger amount of data was collected. Data from the other points 
was corrected based on the data from the "known" point. 

During periods of about eight hours, the receiver was moved from 
one survey point to another. The antenna was set up over the points 
and time difference data was collected for periods of about 20 to 30 
minutes at each point. The time required to move the receiver from one 
point to another varied from about 30 to 60 minutes. 

Data was collected by a Racal Megapulse monitor receiver with 
internal clock accuracy of about 0.02 usec. A Silent 700 data terminal 
was used to collect the data on cassette tapes. Data was recorded at 
intervals of 5 seconds. Later, the data was transferred to 5" disks 
for use on a computer, 

Since there was no way to collect time difference data at two 
points simultaneously, a method of estimating the time differences at 
the known point was developed. This involved using the available data 
to establish linear regression lines of time difference versus time at 
the known point. 

Using the known time of data collection at the other points, an 
estimate of the time differences at the known point could be used to 
calculate a correction for data at the other points at the times of 


data collection. 


A number of computer programs were written to calculate the 
required information for differential corrections. A Correction For 
Three Secondary Stations was developed mathematically and used to 
process the data in an effort to reduce the random error of the time 
difference readings. 

Based on reduction of standard deviations of the data and 
reduced circular error, the three station correction was 
effective to some degree, although absolute error turned 
out to be the limiting factor in the experiment. Reduction 
of absolute error by means of the differential correction 


was significant. 


It. ‘THE LORAN-C NAVIGATION SYSTEM 


A. Basic Description 


Loran-C is an electronic positioning system based on the reception 
of pulsed LF (90-110 kHz) signals from a master transmitting station 
and one or more secondary stations which comprise a Loran-C "chain". A 
brief description of the system will be provided here as a basis for 
later discussions of the accuracy considerations and differential 
methods used in the experiment. 

In the hyperbolic mode, geographic position (GP) is determined from 
time difference (TD) readings provided by a Loran-C receiver. The 
difference in reception time of pulses received from the master station 
and pulses received from each secondary station give the receiver's 
position in terms of a hyperbolic lattice system. Geographic position 
can be calculated from the TD readings using Loran-C lattice tables 
published by the Defense Mapping Agency or computed directly using 
various computer techniques. Nautical charts overprinted with Loran-C 
TD lattices are also available for plotting geographic position 
directly without converting first to latitude and longitude. 

The hyperbolic lattice system is a series of hyperbolas. A pulse 
(actually a group of pulses) is first transmitted by the master 
station. The wave travels to each secondary station of the Loran-Cc 
chain, amd each secondary transmits its own pulse after a time delay 


(called the coding delay) which is unique to each secondary. 


The master pulse is received first by the receiver. A short time 
later the pulses from the secondaries are received. The time dif- 
ference between reception of the master and each secondary is measured 
in the receiver. Each TD identifies one line of position (LOP) in the 
set of LOP's for a master-secondary pair. Figure II-1 shows a typical 
Loran-C master-secondary pair with a set of hyperbolic LOP's. Each LOP 
in this set represents a line of equal time difference between 
reception of the master pulse and reception of the pulse from the 
secondary. Figure II-2 shows an example of a hyperbolic lattice system 
formed by two master-secondary pairs. 

As with latitude and longitude, two LOP's are required to define a 
geographic position. Unlike latitude and longitude, two TD readings 
from different secondaries do not uniquely define a GP, since the TD's 
may intersect at two points within the coverage area. This is normally 
not a problem, since even an approximate knowledge of the geographic 
position of the receiver will allow determination of the correct 
intersection point. 

When more than two LOP's are available to define a GP (when signals 
from more than two secondary stations can be received), the 
intersection of the three LOP's define three geographic points (see 
Figure IV-2). The accuracy of the two best pairs of LOP's is 
sufficient for most navigation purposes, but it is possible to combine 
three or more LOP's to estimate the position of one point. A method 
using corrections from three LOP's is used in the experiment and will 


be discussed in Section IV.B. 


FIGURE II-1. Typical set of hyperbolic LOP's associated with a master- 
secondary pair. 


FIGURE II-2. Typical hyperbolic lattice formed by sets of LOP's 
associated with two master-secondary pairs with a 
common master station. 


B. Time Difference Equations 


The difference between arrival time of the signal from the master 
and arrival time of the signal from a secondary can be described by 
the following equations (Ref. 1): 


The basic time difference equation for Loran-C is: 


lm =-t-t fee 2 | 
where: 


TD 


time difference at a given point P. 


ts, = time of reception of secondary pulse. 


t, = time of reception of master pulse. 


A more useful form of this equation takes into account the 
components which make up the travel time of the signals over each 


section of the paths from the stations to point P: 


T = (Ts + Drs + ED) - (Ty, + DI,,) + ASF [II-2] 
where 
TD = Time Difference at a given point P. 
Ts, = Distance in time units between the secondary station 
and point P. 
T, = Distance in time units between the master station 
and point P. 


DT. = Secondary Phase Correction for all seawater path 
between secondary station and point P. 


DI, = Secondary Phase Correction for all seawater path 
between master station and point P. 


ED = Emission Delay for the secondary station. 


ASF = Additional Secondary Phase Correction for point P. 


Equation II-2 contains two corrections which are applied to Loran-C 
TD readings. The velocity of electromagnetic waves in air under 
standard conditions is changed when the wave propagates over a surface. 
The Secondary Phase Correction Factor (DI, and DI, in equation II-2) is 
a time correction for the change in velocity caused by propagation of 
electromagnetic waves over seawater. This change in travel time can be 


calculated by the following empirically derived equations (Ref. 2): 
For distances greater than 100 statute miles (T > 537 usec): 


DI = Ad/T + Ai + A5*T [ Liss 
where: 


DI = Secondary Phase Correction for all seawater path. 


Ay = 129.04398 psec’ 

A, = -0.40758 psec 

A> = 0.00064576438 

T = Travel time (T, or T, 1n equation II-1) in usec. 


For distances less than 100 statute miles (T < 537 usec): 


DI = Bo/T + B, + BAT [II-4] 
where : 
Bp = 2.7412979 psec’ 
B) = -0.011402 usec 
= Q.00032774624 
T = Travel time in msec. 


In equation II-2 the Additional Secondary Phase Factor (ASF) 
corrects for propagation over land (actually any other non-seawater) 
surfaces, which have a wide variety of conductivities (see Section 
III-C-2 or Ref. 3 for further explanation of ASF Corrections). The ASF 
correction is in addition to the Secondary Phase Correction for all 
seawater path. 

Secondary Phase Corrections for all seawater path are always added 
in before calculation of the geographic position from TD's. ASF 
corrections are much smaller, and are often determined from tables 
which require a dead reckoning position within 10' of latitude and 
longitude of the final position. If no DR position is available, the 
position calculated from TD's corrected only for all seawater path is 
used aS a DR position to determine ASF corrections. When used in 
hyperbolic mode, the ASF correction in Equation II-2 is a net 
correction which applies to the particular master-secondary pair (the 
difference between the range-range corrections for each station). 

Fmission Delay (ED) in Equation II-2 includes a coding delay, which 
is the time between reception of the master pulse at the secondary 
station and transmission of the secondary's pulse, plus one way travel 
time from the master to the secondary, including the Secondary Phase 
Correction for all seawater path for the master-secondary baseline path 
(Ref. 4). 

The result of application of these corrections to the basic time 
difference between reception of master and secondary signals is a 
reduction of the time difference to a path through air with an Index of 


Refraction of 1.000338 (Ref.4:p.4). 


C. Differential Corrections 


Differential techniques can be used to improve the absolute 
accuracy of Loran-C TD readings. Differential Loran-C uses monitor 
receivers located at fixed points in or near the area where improvement 
1s needed. The most common method is to use TD readings of the 
monitor receivers to determine the magnitude of the variations in the 
TD's from the long term average of the TD's at the fixed points. These 
variations are then applied to receiver readings in the area of 
interest before geographic positions are calculated from then. 
Corrections are updated periodically as explained later in this 
section. (Ref.1:p.434) 

The differential correction (including geographic correction) to 
Loran-C readings at an unknown point Py with a monitor receiver at a 


known point P; can be described by the following equation: 


TD at Pp = Received TD at P5 + Differential Cor. + Geographic Cor. 
or 
TD2c = TD2 + (TD, - TDic) + TDg (II-5] 
where 


TDoc = Corrected TD reading at Point P5 (the unknown point). 


TDjc = Long term average TD at P, (location of the monitor). 
TD, = TD reading of receiver at the known point (Pj). 

TD2 = TD reading of the receiver at the unknown point (P92). 
TDg = TD variation which is a function of location. 
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For practical purposes it is possible to use any correction from 
real time to a correction averaged over a period of up to several days. 
The method used in this experiment (described in Section IV-C) was to 
Getermine a regression line of monitor receiver TD readings at the 
known point. A computed correction was applied to the TD reading at 
the unknown point based on the time of the TD reading at the unknown 
point and the equation of the regression line. This method was used 
because only one monitor receiver was available and it was used as both 
a monitor at the known point and at the unknown points. § =In this 
experiment a correction based on reception of three secondary stations 
(rather than two) was applied to the TD readings from each secondary 
station before the regression line parameters were calculated, as 
explained later in Section IV-B. 

The geographic correction - that part of the difference in TD 
readings at points P, and Py which is a due to the differences in ASF 
corrections at the two points - is a particular problen. For 
geographic points at sea, tables of ASF corrections within 10' of 
latitude and longitude are available, and may be accurate enough to 
provide corrections within the necessary accuracy. Tables of 
corrections for geographic points on land are not published, but were 
obtained for this experiment from the Defense Mapping Agency 
Hydrographic/Topographic Center. 

Test results of differential Loran-C on the East Coast U.S. chain 
have shown accuracies of about 4.6 m CEP at distances up to 140 km from 


the monitor site (Ref.1:p.434). 
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TIT. GENERAL ANALYSIS OF ERRORS 


To quantify the amount of error and compare the accuracy of 
different types of corrections made to Loran-C geographic position 
fixes, it is necessary to provide some basic definitions. The 
following sections give the definitions used in this study for the 
types of accuracy and measures of accuracy. One and two dimensional 
sources of error in Loran-C positioning are listed and approximate 
values for the error are given where possible. Error reduction based 


on differential correction is described. 


A. Types of Accuracy 


There are several ways to specify the accuracy of electronic 
navigation systems. The following definitions of three types of 
accuracy have been found useful in this study (from Ref. 5): 

Absolute Accuracy: For purposes of this study, absolute accuracy 
is a measure of the ability to determine geographic position (in terms 
of latitude and longitude or other fixed coordinate system) from a 
given set of Loran-C coordinates (Ref.5:p.172). The error includes all 
types of systematic and random error associated with the system. 

Absolute accuracy in this experiment refers to the accuracy of the 
calculation of the positions of the "unknown" points based on receiver 
TD readings, three station corrections, and differential corrections, 


compared to the known published positions of those points. 
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Repeatable Accuracy: The cammon navigation definition of 
repeatable accuracy is the ability to return to the same point using a 
given set of electronic coordinates (Ref.5:p.172). Error is measured 
in linear units. This differs from absolute accuracy in that the 
position itself need not be known in terms of some other coordinate 
system. Repeatable accuracy is normally better than absolute accuracy 
because the part of systematic error which is a function of position 
has been removed (the position remains the same), leaving only random 
error and time dependant systematic error. Repeatable accuracy in the 
context of this experiment is a measure of the precision of measurement 
of the TD's, which is most affected by the electronic circuitry of the 
transmitters and receivers involved (over short periods of time in 
which terrain and atmospheric factors do not change significantly). 

Repeatable accuracy is Similar to the accuracy expected from 
differential corrections. The difference between the two is that 
Since repeatable accuracy refers to the ability to returm to the same 
point using the same set of received TD readings, there is no 
difference in the ASF (terrain) correction but there may be some 
difference due to variability of TD's with time. 

Differential corrections have different ASF corrections 
at different points. These must be calculated and applied 
to the TD readings. At a single point the ASF corrections 
will probably not change over the time period of interest, 
but over longer periods some error will be introduced. The 
uncertainty in ASF corrections could be a significant 


source of error. 
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Relational Accuracy: The accuracy with which it is possible to 
determine the position of one observer in relation to another observer 
using the same system (Ref.5:p.172). A form of relational accuracy 
enters into differential corrections when two different receivers are 
used. Even at the same point at the same time there will usually bea 
slight difference in the Loran-C readings of different receivers. This 
is caused by errors associated with the receiver circuitry and random 
errors over the slight differences in the time at which the readings 
are made by each receiver, since all other factors are the same. The 
same type of error exists with two different receivers separated by 
some distance (as with differential corrections), although this type of 


error is usually small compared to other sources of error. 


B. Measures of Accuracy 


In addition to the types of accuracy listed above, it is necessary 
to specify how accuracy is to be measured. This can be viewed as an 
estimate of the amount of error. The following are definitions of the 
basic methods used here of specifying the measures of accuracy of 
position. Their usefulness and limitations in this application are 
also described. 

Root Mean Square (RMS) Error: This is the linear (one dimensional) 
error in the location of an LOP (Ref.5:p.172). It can be specified in 
linear units or time units (meters and microseconds are used here). 
The normal specification of RMS error is the number of meters or 
microseconds representing some whole number of standard deviations from 


the mean. In the application used here, the variation in TD readings 
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from each master-secondary pair can be expressed as a standard 
deviation in microseconds, which corresponds to a displacement of the 


LOP by a corresponding number of meters measured normal to that line of 


position. 
Circular Error Probable (CEP): This is the two dimensional error 


which specifies the radius of a circle within which there is a 50% 
probability of being located, based on positions calculated from the TD 
readings of the receiver. The estimated position of the point (based 
on latitude and longitude or other coordinate systems) and the radius 
of the circle are used to specify this measure of accuracy. Circles 
other than the standard 50% circle can be used (67% and 90% are 
common). (Ref.5:p.173) The main advantage of using CEP (rather than 
the error ellipse described below) is that direct comparison between 
accuracies of different systems, techniques, or at different points can 
be made. CEP in this experiment is used as a measure of the precision 
of the TD's indicated by the receiver at a fixed point arn the 
precision of corrected TD's based on these uncorrected TD's. 
Corrections for the absolute position (the bias of the TD's from true 
readings) 1s of secondary importance in this application. 

Circular error estimates in this experiment are based 
on the radius which 50% of the data are greater and 50% are 


less, as calculated by a sorting program. 
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Error Ellipses: Error ellipses are the preferred way of 
describing the accuracy of electronic navigation systems for many 
purposes. As with circular error, accuracy is specified in terms of 
the probability of an estimate of geographic position falling within a 
geometric figure centered on the actual geographic position (an ellipse 
rather than a circle is used) (Ref.6). Unlike the circle used in CEP, 
an ellipse provides additional information about the directional 
Characteristics of the accuracy, which can be very important in som 
applications. The specifications for elliptical accuracy are usually 
the lengths of the major and minor axes and the orientation of the 
ellipse (such as the azimuth of the direction of the major axis). A 
disadvantage of using error ellipses in this experiment is the 
difficulty of directly comparing the accuracy of different ellipses 
with different orientations. For that reason most of the accuracy 
comparisons in this study have been specified in terms of circular 
error. Error ellipses have been used where appropriate. 

Error ellipses are related to the orientation of two intersecting 


LOP's by the following equations (Ref.6:p.72-74): 


04° sin 2p 
Can: 26.03 eee (Tite) 
01% (cos 2p) + 052 


042 ata 05° + ((042 + G5°)¢ = 4(sin2$) 012052) % 
Oye = — sss pec 
2 sin-¢ 
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. l. 
012 + a5% = ((012 + 09%)% - 4(sin*¢) 01202) # 


PN, EAS = = 
Oy =a (ILI-3] 
2 sin¢¢ 
where: 
@ = the angle of the major axis of the ellipse 


from the direction of the LOP with smaller 
variance, toward the LOP with greater variance, 
in the direction of the smaller angle. 


0)“ = variance of the LOP with smaller variance. 
05" = variance of the LOP with the greater variance. 
® = the smaller angle between LOP's. 


These equations were used to calculate the dimensions and 
orientation of the error ellipse for Range-7, the "known" point used 


for differential corrections (see Figure V-3). 


C. Sources of Linear Error in Iorarmc Tine Differences 


The position fixing accuracy of the LORAN-C system depends on a 
number of factors. These can be divided into two categories 
(Ref.7:p.18): temporal errors affecting the repeatability of the 
system (errors which vary with time), and systematic (time invariant) 
errors. Both of these types of error may also be functions of 
“geographic position. 

The effect of these errors is to produce TD readings at the 


receiver which differ from those which would be expected on a perfect 
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hyperbolic lattice on the surface of a perfect ellipsoid, resulting in 
error in the calculation of the geographic position. A linear error of 
1.0 usec in time difference is approximately equivalent to 150 m on the 
baseline, with one dimensional error at other points being up to about 
ten times greater, depending on the location within the geometry of the 
system. 

The main sources of error in the determination of TD's used to 
calculate geographic position can be classified as follows (Ref.8): 
Transmitter synchronization timing errors. 


Errors associated with propagation effects. 
Atmospheric noise. 


Receiver errors. 


1. Transmitter Synchronization Timing Errors 


Since time differences between reception of a master pulse and 
secondary pulses are used to calculate geographic position in 
hyperbolic Loran-C, it is essential that the synchronization of timing 
between transmissions from the master and secondaries be highly 
precise. 

In order to accomplish this, each station is individually timed 
using cesium clocks which drive pulse amd group repetition interval 
timing circuits. Independent monitor receivers are used to provide 
data to computers which adjust the relative timing between master and 
secondaries of the chain in small steps (usually 10 to 20 ns). 
Absolute timing of the master stations is adjusted to operate in 


coordination with Universal Time, Coordinated. (Ref.8:p.1129) 
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Estimates of timing error due to transmitter synchronization 
are on the order of 0.05 usec (Ref.9:p.537). Actual data from an 
operational Loran-C chain showed time differences controlled within 
0.024 wsec of the standard value, with a standard deviation (in 
absolute time) of 0.032 psec (Ref.10:p.227). A transmitter 
synchronization error of 0.01 usec is believed to be attainable if 
ground monitor data is used to remove the mean synchronization error 


(Ref.11:p.401). 


2. Errors Associated with Propagation Effects 


This includes most types of error external to the electronics 
of the system, and is the main source of error in the determination of 
geographic position. 

Factors which influence groundwave propagation include the 
following (Ref.12:pp.39-53 and Ref.13:pp.173-187): 

Electrical Conductivity of the surface material 
Dielectric Constant of the surface 


Index of Refraction of the air 
Lapse Rate of the Index of Refraction with altitude 


These factors are complex functions of several variables, the 


main ones being: 


¢ type and composition of the surface over which the wave travels 


e topography of the surface 
¢ meteorological conditions 
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At any particular geographic position the effect of these 
variables may or may not be considered a function of time, depending on 
the time periods of interest. 

Also of importance is the effect of these factors on TD as a 
function of geographic position. For each different geographic 
position, the electromagnetic waves from each transmitter must traverse 
a different route, with different types of surface, topography, and 
atmospheric conditions, to reach the receiver. 

The effects of irregular terrain are practically time 
invariant, while ground conductivity and dielectric constant may vary 
seasonally (due to changes in moisture content of the soil, changes in 
vegetation, or other factors). Observed cyclic annual variations on 
the order of 0.5 wsec are probably due to seasonal variations in ground 
conductivity and dielectric constant and to seasonally varying 
atmospheric conditions. 

Meteorological conditions sufficient to cause noticeable 
variations in TD readings may change in very short time periods (days 
or hours). Short term variations due to meteorological effects (cold 
front and warm front passage through the line of propagation) on the 
order of 0.025 usec have been observed. These correlate with changes 
in the dry adiabatic lapse rate of the index of refraction. 
(Ref.14:p.1119). 

For most differential Loran-C purposes it is necessary to 
consider only those temporal propagation effects which occur in short 
time periods (hours or days at most, for most applications). The 


variation of propagation effects with geographic position must also be 
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considered for precise positioning. Many of these variations can be 
corrected by a computation of the Additional Secondary Phase Factor 
(ASF) . 

The Additional Secondary Phase Factor (ASF) is a measure of the 
systematic error resulting from differences in velocity of electromag- 
netic waves over ground with different conductivities and dielectric 
constants. The Defense Mapping agency uses a computerized system of 
maps with estimates of ground conductivity (effectively including 
dielectric constant) over large areas, amd calculates ASF corrections 
in the Loran-C coverage areas. These corrections are published as 
tables of ASF corrections to be used in correcting Loran-C TD readings 
from receivers. 

ASF corrections can be as large as 4 microseconds (Ref.3). 
The corrections provided in the Defense Mapping Agency tables require 
estimates of the geographic location of the receiver to the nearest 10' 
of latitude and longitude. ASF corrections based on the known 
positions of the survey points used in this experiment were provided by 
DMA using a computer program which calculates the corrections using the 


Millington Method (Ref.3:p.V). 


3. Atmospheric Noise 


The effects of atmospheric electrical noise are reduction of 
the range at which Loran-C signals can be effectively received and an 
increase the error in TD's computed by receivers in marginal reception 


areas. This is primarily dependent on the circuitry of the receiver. 
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In marginal areas the effect of atmospheric noise will be an 
increase in the standard deviation of the TD's and increased difficulty 
in locking on am maintaining reception of the weaker signals. 
Atmospheric noise varies with time and geographic location, with a 
Significant diurnal variation (noise levels are higher at night). In 
this experiment, significant effects due to local sources of electrical 


noise were noticed in the Monterey harbor area. 


4. Receiver Errors 
Receiver errors are of two main types: 


¢ timing error caused by limitations of the internal clock 
Circuitry of the receiver 
* errors caused by processing time of the receiver 
For a receiver at a fixed geographic location, processing time 
is not a source of error. The receiver will be at the same location 
regardless of the time required to output the TD's or geographic 
coordinates. The main source of receiver error at a fixed location 
(other than the noise considerations discussed above) is the timing 
error in the internal clock circuitry. lLoran-C receivers vary greatly 
in their timing errors. The Racal Megapulse Monitor receiver used in 
this experiment has a specified timing error of 0.02 psec or better. 
For mobile receivers the processing time of the receiver 
circuitry can be important. After the signal is received it must be 
processed and either stored in the form of electronic information or 
visually displayed. Error is introduced if the receiver moves during 
this time interval. 
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Obviously, the amount of error in the determination of 
geographic position from this source depends on the processing time and 
the speed of the receiver. Receivers which output data in the form of 
geographic coordinates mst calculate latitude and longitude from the 
TD readings of the receiver, requiring more processing time an 
resulting in greater error from this source. 

At the speeds normally used in hydrographic” surveying 
(typically less than 20 kt) the processing time is insignificant for 
practically all receivers. For data which is stored in the form of 
receiver output TD's for later correction and conversion, as with the 


Racal Megapulse receiver used here, error should be insignificant even 


at much higher speeds. 
D. Geometric Dilution of Precisign 


This source of error in geographic position results from the 
crossing angle of the Loran-C LOP's (lines of position). Figure Iv-1 
shows two hyperbolic LOP's with standard deviations 0, and oy (in time 
difference) and crossing angle a. The optimum crossing angle of any 
two LOP's is 90°, resulting in minimm error in geographic position for 


a given set of standard deviations. 
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E. Summation of Errors 


1. Absolute Accuracy vs. Precision 


The absolute accuracy of an electronic positioning system can 
be viewed as a measure of its ability to give the geographic position 
of the receiver. Precision is a measure of the variation of individual 
time difference readings. If the receiver is set up at a point fora 
sufficient amount of time, many of the random errors can be averaged 
out; the average of a large number of readings can provide a basis for 
determining absolute accuracy. 

In practice this is more difficult than it would seem, since 
some factors which do not vary significantly over short periods of time 
(hours or days) may vary Significantly over longer periods (weeks or 
months). It is necessary to consider the time periods involved in 
order to determine which factors to consider random and which to 
consider "bias" in the system. 

For this experiment, the sources of random error are considered 
to be transmitter timing errors, some meteorological conditions, 
atmospheric noise, and most receiver errors. Errors caused by long 
term "bias" of the system are due primarily to ground conductivity and 
type of surface, some meteorological effects, and any non-random bias 
in the timing of the transmitters. 

Taking these factors into consideration, data taken at the 
"known" point (Range-7) of the differential system was processed to 
remove random error by averaging over a period of time, while long term 


error waS removed by comparison of the TD's to their expected values 
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derived from calculations provided by computer computations of the TD 
at that point. 

Errors which may change within the time frame of the experiment 
were reduced by a linear regression of the trend of data taken at the 
known point and projected to the correction at the unknown point 


mathematically. 
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IV. THREE STATION AND DIFFERENTIAL CORRECTIONS 


A. Basic Geometry of LoranC Hyperbolas 


This section provides a summary of the information necessary for 
the analysis of Loran-C corrections for three secondary stations and 
for differential Loran-C corrections presented in parts IV.B ard IV.C, 


adapted from Laurilla (Ref.1) and others as noted. 


1. lLanewidth 


The baseline of a Loran-C master-secondary pair is the shortest 
path over which the transmitted signals may travel between the master 
station and the secondary station. On a spherical earth model it is 
the arc of the great circle through the positions of the master and 
secondary stations. On an ellipsoid it is the geodesic from the 
master to the secondary station. In the set of hyperbolas of time 
difference for the master-secondary pair, the separation between 
hyperbolas is a minimum constant value along the baseline. 

A lane is the distance between two adjacent hyperbolas of unit 
time difference. The unit normally used is the psec (1076 sec). Ina 
series of hyperbolas, as shown in Figure II-1, it is obvious that the 
lanewidth varies with location. At any point,  lanewidth can be 
viewed as the inverse of the gradient of the ratio of the change in 


time difference to the change in distance (Ref.15), or: 


es | Cetemysit | [rvaig 
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Figure IV-1. Typical intersection of two lines of 
position. Most probable point is at the interséction of 
the two LOP's. 


LOP} 
LOP> 
LOP3 
maqure -IV-2. Typical intersection of three lines of 
position. When three or more secondary stations can be 


received at a geographic position, calculated position 
depends on the LOP's which are used to make the 
calculation. 
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The lanewidth at any point in the set of hyperbolas can be 
expressed in terms of the lane width on the baseline. Baseline lane 


width can be expressed by the following equation: 


Cc 


ue [Iv~2} 
where: 
Ix, = Lanewidth on the baseline (meters/pusec) 
c = Speed of electromagnetic waves on the baseline (m/usec) 


(note that time units are in the unit lanewidth) 


Lanewidth at any other point (not on the baseline) can be 
expressed by equation IV-3a, where #8 is the angle subtended by the 
baseline from the point. Note that this angle is the absolute 
difference between the azimuths from the point to the master and from 


the point to the secondary station. 


L = Ip{sin(@/2))71 [Iv-3a] 
or 
Iia= she 
where: 
L = lLanewidth at a point (per unit of time) 


Iy, = Lanewidth on the baseline (per unit of time) 
6 = Angle subtended by the baseline from the point 
G 


= The Lane Expansion Factor, G = {sin(8/2))71 


Note that the lane expansion factor depends only on the angular 


separation of the master and the secondary stations. This obviously 
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depends on the distance of the point from the baseline center and the 
angular separation of the baseline and the vector from the baseline 
center to the point. 

Equation IV-3a provides a means of calculating the normal 
distance between any two hyperbolas, provided they are near enough that 


their lanewidths are not significantly different: 


d=It [IV-3b] 


where: 


d = distance in meters between the two hyperbolas. 


ee 
ll 


lanewidth at the point of interest (in m/psec). 


Gr: 
li 


time difference between the two hyperbolas in usec. 


2. Direction of Hyperbola 


The direction of a hyperbola at any point on the hyperbola is 
along the tangent at that point, directed away from the baseline of the 
master and secondary stations used to define the hyperbola. For 
applications in this experiment, Laurilla (Ref 1) provides the follow- 
ing useful definition of the direction of a hyperbola: " .. the 
direction of a hyperbola at any point coincides with the bisector of 
the angle formed by the lines joining the point to the pair of 
stations." (Ref.1:p.92) 

For purposes of this experiment, it is convenient to express 
the direction of the hyperbola by its azimuth. This azimuth can be 
expressed in terms of the azimuths (from north) of the master and 
secondary stations by: 
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(Qi, + As) /2 (where |a,,-a.|>180°) [IV-4a] 


aH 
li 


180° + (Q_ + @,)/2 (where |a,-a,|<180°)  [IV-4b] 


=i 
iH 


where ; 
@, = azimuth of the direction of the hyperbola. 


@, and a, = azimuths of the master and secondary from the point. 


3. Applications of Geometry of Hyperbolas to Error Analysis 


The geometric description of Loran-C hyperbolas in parts IV.A.1 
and IV.A.2 is particularly important in the analysis of error and the 
derivation of corrections used in this experiment. Some of the direct 
applications are outlined in the following paragraphs. 

The error expressed as an error ellipse requires the directions 
of hyperbolas in order to determine the orientation of the ellipse. 
The dimensions of the ellipse require the conversion of variances of 
time differences of each hyperbola to variances of distance in the 
major and minor axis of the ellipse. Standard deviations of time 


differences can be converted to distances by the following equations: 
0g = O¢L { TV—-5a } 


Gg’ = 04°L? [IV-5b] 
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The correction for three secondary stations uses variances in 
time differences converted to distances as an input, using equation 
IV-5b. Converting time difference to distance in this way correctly 
weights the precision of the LOP's when determining the correction 


factor based on geometric considerations. 


B. Correction for Three Secondary Stations 


In situations where three secondary stations can be received, the 
result is three LOP's which do not, in general, intersect at one point 
(see Figure IV-2) when they are used to calculate the geographic 
position, although they must in fact intersect at the common point at 
which the readings were taken concurrently. Under usual navigation 
conditions the two best LOP's are used to determine a geographic fix. 
A drawback of this method is that it fails to consider all of the 
available information (i.e. one of the available lines of position is 
not used in the calculation of the geographic position). By utilizing 
all three LOP's it should be possible to improve the accuracy of the 
geographic fix. 

Sections B.1 amd B.2 will develop a mathematical adjustment which 
can be used to calculate the most probable corrected TD readings of a 
Ioran-C receiver used as the fixed receiver in a differential Loran-C 
system utilizing three LOP's (i.e. receiving the master and three 
secondary stations). When the correction is applied, the three 
corrected TD's will intersect at the most probable point, based on long 
term average values of the TD's. The derivation of this adjustment 


could take any of several possible forms, but a primary consideration 
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here is the need for a form which is easily applicable to computer 
processing. Section B.3 describes a method of applying the corrections 
for non-differential use, using only one receiver and with no known 
fixed point. 

Section B.4 will develop regression equations for calculating the 
corrected TD's directly from time for use as a differential correction 
for a remote receiver. Note that the methods developed in these 
sections are primarily for use in the differential Loran-C computations 


necessary for this experiment. 
1. Geometric Considerations Used to Develop a Correction Factor 


The following assumptions are made for this derivation: 
e Each LOP has a unique direction (described in Section IV.A.2) 


¢ Long term average TD's (designated TDl,,, TD2s,, and TD3,,) for 
for the known point are available 


e Variances of the TD values (designated Or, On"; and 03°) are 
known or can be calculated 


¢ Variations of TD's from their average values is in the form of a 
Gaussian distribution 


« Over the relatively small area being considered, the surface can 
be considered a plane and LOP's for the same master-secondary 
pair can be considered parallel straight lines 

Several things should be noted about these assumptions. Long 
term averaging times for this experiment will be on the order of a few 
hours during acquisition of data; for purposes of the derivation, 
standard deviations are converted between time units (usec) and 


distance units (m) using Equation IV-5a; and, the assumptions of a 
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FIGURE IV-3. Geometry used in the derivation of a correction factor 
for three LOP's not intersecting at the same point. 
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Gaussian distribution, plane surface, and parallel LOP's are reasonable 
for this type of data (Ref.6:p.8). When two IOP's with different 
standard deviations intersect (see Figure IV-1), the point of 
intersection is the most probable geographic position (Ref.6:p.8-12). 
The probability of the location of the true position is best described 
using an error ellipse as described in Section III-B. 

A typical situation is shown in Figure 3, a set of three 
concurrent LOP's intersecting at three different points. Point P is 
the location of the monitor receiver used to determine differential 
corrections, with TDlsy,, TD2s,, and TD3,, being the long term average 
TD values at that known point (note that in figure IV-3, the LOP's of 
TDlay, TD2ay, and TD3,, are not shown, but are parallel to TD), TD, 
and TD3, respectively, and pass through point P, as shown in Figure 
IV-4). Although the long term average LOP values intersect at point 
P, triangle P)9P}3Po3 represents a triangle of fixes formed by the set 
of TD values which may occur from any concurrent set of TD's, which are 
unlikely to equal the long term average. Each set of TD readings 
(which occurs in the data every few seconds) will result in a new 
triangle of fixes. 

The vectors rj, Yo, and r3 represent the difference between the 
individual TD values in a concurrent set of data and their respective 
long term TD values. The r vectors are at right angles to both the 
long term and current LOP's. For example, if the long term average of 
TD, (TDlay) is 29939.200 and the value of TD, in a set of TD readings 
currently being considered 1s 29939.300, the value of rj for the 


current set of TD readings is 0.100 (all in usec). When converted to 
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meters using equation IV-3b, this is the distance between the average 
TD values (LOPiz,) and current TD values (TDj|). These LOP's can be 
considered as parallel straight lines over the relatively small 
distances considered for this correction. Note that there are three 
LOP's in each set of concurrent data representing TD values for three 
master-secondary pairs. The master station is common to each of these 
pairs. 

Within the triangle of fixes P}2P,3P53 formed by the three 
LOP's of a set of data, point P' represents the most probable point of 
the true geographic position indicated by that particular concurrent 
set of data compared to the long term average data. The distances of 
P' from each IOP (the lengths rj') are related to the relative 
precision of each LOP, P' being closer to LOP's with higher precision 
(to be discussed in more detail in the next section). 

The problem must be solved in general terms so that any set of 
TD's in the data can be corrected to a new set of TD's which intersect 
at the most probable point. Since the triangles formed by all sets of 
data are similar triangles, the lengths of r,', ro', and r3', will have 
the same proportion for each set. The corresponding rj' values between 
sets of data will be related by a common factor, which will be 
designated as the Correction Factor, k. 

Other information which should be noted about Figure IV-3 
include the following: 

¢ The following sets of lines can be considered parallel over the 
area represented by the figure: TDi,y || TD] ; TD2ay || TD; 
EDI has Fy |. r= XD x5 and 13 i na. 
(continued on next page) 
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* The following sets of lines are perpendicular: r, andr ,' 1 
TD1,y and TD, ; r> and r5' 1 TD2,, and TD, ; . r3 andra 
TD3ay and TD3. 

¢ The azimuths (from north) of the directions of the IOP's 
represented by TDlz,, TD2ay, and TD33, are aj, a2, amd az, 
respectively. These are the same as the azimuths of TD), TD>, 
and TD3, respectively. 

* The following angles of triangle P}2P)3P23 are defined by: 

Q19 = ap - ay 3 a3 = 1/2 — (a3 — a) ; and a93 = 23 —- Gy. 
These angles will be the same in the similar triangles formed 
by any lengths of rj, ro, and r3. 

- IOP's 1 and 3 have the largest absolute angle (+ or -) between 
them, and LOP 2 is the intermediate angle. 

The next section will provide a mathematical derivation of the 
expression for a correction factor to determine the lengths of rj', 
ro', and r3' and the corresponding time adjustments to TD), TD), 


TD3, based on the geometry presented in this section. 
2. Mathematical Derivation of the Correction Factor 


An important question in this analysis is how to weight the 
relative precision of each LOP in such a way that it can be used in 
the calculation of a correction factor. Reference 16 shows that for 
this type of data the "weights are inversely proportional to variance" 
(Ref.16:p.60), meaning that the lengths of the rj;' corrections should 
be proportional to the variances (in units of length) of : the 
corresponding TD values for each LOP. Using units of length, in the 


general solution for all sets of data, the following equation applies: 


ry'Ly = kop*Li° [IV-6] 
where: 
rj' = corrections for TD's 1, 2, amd 3 (in psec). 
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a:2 = variances of TD's 1, 2, and 3 (in usec?) . 


nw 
li 


the Correction Factor (units are meters) 
L; = lanewidths of TD's 1, 2, and 3 (in m/usec). 


Note that 0;°L:- is a unitless weighting factor. 


For any triangle there will be some unique value of k which 
relates the distances from point P' to each LOP to the variances of the 
LOP's according to equations IV-6. To abtain a general expression for 
k which can be easily calculated by computer, it will be helpful to use 
only the known values of the TD's, variances, azimuths, and lanewidths 
for each LOP. The r values are readily obtained from each set of TD 
data by subtracting the long term average values from current TD 
values. Values of a. and L; are constants for this small area. 

If each LOP is considered a straight line and a cartesian 
coordinate system is assumed with north as the +y and east as the +x 
direction of the coordinate axes, then each LOP can be specified by the 


general equation for a straight line: 
y = mxt bj (gaa, 2,.0r 3) [IV-7] 


Three such equations can be written (one for each LOP), giving 
three linear equations in two unknowns, x and y. These equations will 
not have a common solution when the three LOP's do not intersect ata 
common point. 

For each master-secondary pair, all LOP's of the set will have 
the same slope (m;) within the small area being considered (TD; is 
parallel to TDi,,, for any TD;). The y-intercepts (the b;'s) will be 
zero for the average LOP's intersecting at point P. Y-intercepts for 


oy 


other LOP's with non-zero r values can be calculated from the 
relationships between azimuths (denoted aj) and the r values of each 


LOP by the following equation: 


oa les 
4S = ——— (i=1, 2, or 3) [1v=8] 
cos (a@;-71/2) sin aj 
N (y) 


Set ee te ee ee 





E (xX) 


FIGURE IV-4. Typical average LOP through the origin of a cartesian 
coordinate system and corresponding parallel LOP of an individual 
TD reading. Slopes are the same for both. Y intercept is the b 
value in equation for the LOP (zero for the average LOP). 
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Slopes of the TD's can be expressed in terms of ¢ as: 
m; = -tan ¢; [IV-8b] 


Figure IV-4 shows the assumed coordinate system along with the 
LOP's for which the slopes and intercepts are described. When 


corrections are applied, the corrected r values are as follows: 
a Tee ko: 2 .2 Iv-9 
richi ribi poli [ } 


Note that in the equation for corrections (IV-9) the units of 
Yjcl; are meters. If o;°Li? is considered a unitless weighting factor, 
then the units of k must be meters. This agrees with the units of 
Equation 13. 


The value of r; 1s calculated by: 


rj; = TD; - TDiay [Iv-10} 
where 
rj = deviation of an LOP from the long term average (usec). 
TD; = time difference reading of a master-secondary pair at 


pone <P. 


TDigy = long pee of the TD for the same pair at 
point P. 


Figure IV-3 shows that the rj;' corrections must be subtracted 
from r, and r3, amd added to r5, for the configuration shown. It can 
be demonstrated that the r,' corrections are always of the same sign 
(both + or both -) for the two most widely spaced LOP's and of the 
opposite sign for the intermediate LOP, where the direction vectors at 
point P are in the positive directions of increase of TD and the angle 
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is in absolute terms (maximum being 180°). This is true for all 
positive and negative r values. Positive and negative k values will 
reverse the signs of all corrections in equations IV-9 so that they 
will be correct for all configurations. 

The corrected r values of equations IV-9 can be substituted 
into equations IV-8, providing a set of y-intercepts (the b; values in 
the set of general equations IV-7) for the corrected LOP's which 
intersect at point P'. These intercepts contain the unknown Correction 
Factor k. Note that the slopes of the corrected LOP's are equal to 
those of the corresponding uncorrected LOP's, meaning that the mj; 
values in equations IV-7 are unchanged. 

Substituting the corrected y-intercept (bj) values into 
equations IV-7 results in three linear equations in three unknowns: x, 


y, amd k: 


rly - ko,*L,? 


Yo = (tan hy) ee eee (IV-11a] 
COS $7 
rgly + ko2*Lp* 
Y =" (tan ©). .. eee [IV=iiiay 
COS $5 
r4L3 ~ ko37L,¢ 
Y= (Can ds) x + ay [IV-11c] 
COS $3 
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Equations IV-11 can be rewritten as the following set of linear 


equation in the unknowns x, y, am k: 














04°I4° rly 

(tan ¢))x + y + ( i — [IV-12a] 
02°*Ly* r2lo 

(tan ¢do)x + y - ( )k = —_——— [IV-12b] 
03°1,° r3l4 

(tan ¢3)x + y + ( S$ [IV-12c] 
COS $3 COS $3 


Solving equations IV-12 by Gaussian elimination gives the following 


general expression for the Correction Factor k: 


tan $3 - tan ¢> tan $3 - tan a) Bap @oe = eam) 
J | hae. 1-22) om oo 335 °° 
: COS $) COS p2 COS $3 
tan $3 - tan > Ay yee) (7. Came = can 7 
0,2 “+ 205°  ___—_—— + L, “0 ,>@ ______—_ 
COS $1 COS $2 COS $3 
[IV-13] 
where: 
r; = differences from average TD's (}isec). 
k = Correction Factor for three master secondary pairs (m). 
gj = directions of increasing TD (azimuth from north). 


oa = variances of LOP's (part of unitless weighting factor). 
L; = lanewidths (in meters/p"sec) 


Note that ihe eh terms in denominator are unitless weighting factors. 
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Equation IV-13 may be simplified by the following substitutions: 


Ay = Iy(tan ¢3 - tan ¢2)/(Cos 4) [IV-14a] 
Ap = Lg(tan $3 - tan ¢))/(cos ¢9) [IV-14b] 
Ag = La(tan ¢2 - tan ¢1)/(cos $3) [IV-14c] 
Ay = 19017A, + Ip097A> + Ly037A3 [IV-14d] 


Resulting in the following equation for k: 


k = (Ajry — Agro + A3r3)/Ag [IV=159 


The x and y coordinates of point P' can be calculated by: 


Y1COS $> - roCOS ¢) - ko, “cos a ko52cos $2 
EEE eee [IV-16] 
(COS $1) (COS >) (tan 2 - tan ¢}) 


- (tan a1)x + (ry, - 017k)/(cos ¢7) [IV-17] 


= 
lI 


Two other useful parameters can be calculated from the values 
of x and y. These are the azimuth and range of the of point P' from 


point P. These are given by the following: 


Azimuth (from north): 
Qppr = 7/2 - tan l(y/x) (x > 0) [Iv-18a] 
apP: 
Note: when x=0 the azimuth is zero for positive y and 7 for negative y 


3r/2 - tan -l(y/x) (x < 0) [IV-18b] 


Range(P to P'): Rpp: = (x2 + y2)% [IV-19] 
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Using the computed value of k for a set of concurrent TD data, 
a correction can be made to the TD values using equations IV-9. Since 
Iv-9 is in meters it must be converted to usec using equation IV-3b. 
The resulting corrected TD's will intersect at the point P’. 

Equations IV-15, IV-16, and IV-17 do meet the requirements 
stated at the beginning of this section. The inputs are limited to 
azimuths, variances, long term average TD values, amd concurrent TD 
data, and the values of k, x, y, and the corrected TD values can be 
computed directly from these by equations which can be calculated 


easily by a computer. 


3. Non-Differential Correction for Three Secondary Stations 


When using the method of Section IV.B.2 with a single Loran-C 
receiver (i.e. not using differential methods) a correction can be 
calculated based on existing computer programs. While not used in this 
experiment, the method of correction will be outlined in this section. 

Since variation from a known point cannot be used, the two best 
LOP's can be used as the basis for this calculation. While not as good 
aS uSing a known fixed survey point, this will provide a correction 


based on three IOP's. The method is outlined below: 


1. Using the two LOP's with the best fix accuracy (available from 
existing computer programs), calculate the geographic 
position of the point of intersection. This corresponds to 
point P. 


2. From the geographic position of P, calculate the TD of the 
third master-secondary pair (the one not used in the previous 
step). 


(continued on next page) 
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3. Using point P calculate the azimuths (from north) of the LOP's. 
Number the TD's as described in the previous sections, with 
the intermediate TD as TD) and the TD's with the greatest 
angular separation as TD, and TD,. 

4. Use TD variances calculated over short time periods in the 
general area, or use the values of metersS/microsecond cal- 
culated using equation IV-3 as the approximate weights of 
each LOP. Substitute these for the o“ values in equation 
LV=13? 


5. Assume r values of zero for the LOP's with the best fix 
accuracy used in step one. 


6. Calculate the difference between the observed third TD ard the 
calculated value obtained in step 2. Use this for the third 
r value. 
7. Calculate corrected r values as in the previous sections. 
8. Using any two corrected TD's, calculate geographic position. 
This is the position corrected for three secondary stations. 
Note that all of the above calculations are based on computer 


methods and programs already in use (Ref. 2) allowing corrected posit- 


ions to be calculated readily using the method outlined. 
C. Differential Corrections 


One of the main objectives of this experiment is to apply 
differential corrections obtained at a known point to time difference 
data taken at unknown points. One problem encountered was the 
availability of only one Loran-C monitor receiver, making simultaneous 
data collection at two points impossible. Differential Loran-C does 
not require continuous updating of the differential correction for 
general navigation purposes (hourly or daily corrections are possible), 
but the need for the most accurate corrections possible for the 


experiment required a continuous update of the correction. 
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1. Mathematical Description of the Differential Correction 


Data from the experiment consisted of a series of readings at 
each point. Data was recorded at intervals of 5.0 seconds for mich of 
the experiment, although other time intervals were also used. Data 
consisted of time versus TD for up to four possible master-secondary 


pairs and the signal to noise ratios for each, in the following format: 


Time (GMT) GRI 1D, TD) 1D; TT, S/N SW Sy Sy 


Only three master-secondary pairs were receivable in the area 
of the experiment. One survey point (Range 7) was designated the 
"known" point, from which differential corrections were calculated for 
data taken at the other points. Over a period of several days, data 
was collected at the known point for periods of 20 to 30 minutes at 
irregular intervals spaced at about one to three hours (data was 
collected at the other points in the intervening periods as described 
in Section A.2). 

To provide estimates of the continuous TD readings at the fixed 
point during all time periods (including the periods between data 
collection at that point), a least squares curve fitting process is 
used. The result is three equations: TD as a function of time, for 
each master-secondary pair. 

These three equations can be used to provide an estimate of the 
TD values at the fixed point at any time. The difference between the 
calculated TD at a certain time and the long term average TD is the 
Differential Correction at that time. The ASF corrections at the fixed 
point and an estimate of the ASF correction at the unknown point (as 
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described in Section II-B) also enter into the Differential Correction 


at the unknown point. An equation is written for each LOP: 


Cy = TDo - TDay + (ASF, - ASF) [IV-20] 
where: 
Cqg = Differential Correction at the unknown point. 
TD. = TD (at known point) calculated from regression equation. 


TD3y = Long term average TD at the known point. 


ASF = ASF correction at the known and unknown points. 
2. Least Squares Regression Equations for Time Difference 


Linear regression equations for each TD can be written in the 


following form (Ref.17): 
TD = Ho sts byt [ FVe2a5) 


The b; coefficients of the regression equation can be found by 


solving the following set of linear equations: 


nbo + by2tj = 2TD; 
[IV-22a] 
Bott; + by=t;2 = 2Xtj;TD; 


An alternate method of calculating the b; values directly is: 


n=tj;TD} - (2tj) (2TD}) 








ae 2 2 
mei 7 ad (IV-22b} 
2TD} 2TD} 
Bo = =a7 
nN nN 
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In equations IV-22, <Xtj, Xt;TD;, etc., refers to the sum of 
the time values, the sum of the products of time and time difference 
values, etc., for the whole set of data available at the fixed point. 


The total number of sets of TD vs t data is n. 
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V. RESULTS AND CONCLUSIONS 


Part A of this chapter presents the results of calculations based on 
the theoretical principles previously discussed and data obtained from the 
experiment. Part B provides a summary of the results obtained, how these 
results differ from what would be expected based on the theory, possible 


reasons for any variations from expected outcome, and a general conclusion. 


A. Calculations 
1. Azimuths, Directions of Hyperbolas, and Directions of TD's 


Azimuths of the Loran-C stations from the survey points are one 
of the basic requirements of several calculations used in this 
experiment. Station positions of Range-7 and Luces Point (the primary 
stations used in this experiment) and other points were obtained from 
Horizontal Control Data published by the U.S. Dept. of Commerce, Coast 
and Geodetic Survey (see Appendix A). Station positions of the Loran-C 
transmitters were obtained from Loran-C Data Sheets provided by the 
Defense Mapping Agency, Hydrographic/Topographic Center (see Appendix 
B). For convenience in avoiding datum conversion, the North American 
Datum 1927 (NAD 27) was used because published positions based on this 
datum were available for all of the station positions used in the 


experiment. Station positions are summarized in Table V-1. 
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Azimuths: Azimuths from Range-7 amd Luces Point to Loran-C 
stations 9940-M, W, xX, and Y were calculated using a Fortran 77 
subroutine provided by DMAHTC (Appendix C). Operation of the. 
program was checked using known data points and azimuths. Results of 
these computations are presented in tabular form in Table V-2. For 
convenience in visualizing the angular relation- ships, results are 
also presented in graphic form in Figures V-la and V-2a. 

Directions of Hyperbolas: The directions of hyperbolas for 
9940-W, X, and Y were calculated using equations IV-4 and the azimuths 
obtained as described in the previous paragraph. A sample calculation 


follows: 
Hyperbola Direction for 9940-X at Range-7: 


(Azimuth of 9940-X) + (Azimuth of 9940-M) 
KDpy = 
Pe 
343.071686°+ 38.358552 
a ee ee ee «=— LOR TT LOLS 


2 
(from North) 


Decimal degrees are used in the calculations for convenience. 
Results of the calculations are again presented in Figures V-1b amd 2b. 

Direction of Increase of TD: Directions of increasing TD are 
at right angles to the hyperbola directions amd increase in the 
direction of hyperbolas which approach the master station. This 
involves either adding or subtracting 90° from the direction of the 
hyperbolas. For checking purposes, a Loran-C chart of the area was 
used. Results are again presented in Figures V-1 amd 2. 
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TABLE V-1. Station Positions (NAD 27). 


Range 7 36° 39' 02.47787" N | 121” 49° 08.58202™ 54 
Iuces Point 36°938' 102524. 121° 55° 38.399" 
9940-—M 39° 33' 07.046" N 118° 49' 52.241" W 
9940-W 47° 03' 48.594" N 119° 49' 52.241" W 
9940-X 38° 46' 57.472" N 122°*..29' 40.0507 
9940-Y 35° 019" 18.34200N 114° 48' 13.946" W 





TABLE V-2. Azimuths From Range-7 and Luces Point. 


9940-—M 21' 30.784" 38 .358551° 
9940-W 920922 5625. 7.53967 38 
9940-X 04' 18.070" 343 .071686° 


9940-Y AT 3G nee on 101.7933975 


08.262" 38.852295° 


NS Gi ye Oe ly oe 


49.981" 344.397217° 


6. oios 101.604533° 





50 


S940-W 72539673 


9940-X 343.071686° 
9940-—M 38.358551° 


SS40—7 ~ 101 793397 = 


FIGURE V-la. Azimuths From Range-7 to 9940-M,W.X,and Y. 


9940-Y 250.075974° ~~ “A 


9940-W 202. 949112" / /| 
9940-X 190.715118° 


FIGURE V-1b. Directions of Hyperbolas at Range-7. 


9940-Y 340.075974°* 


9940-X 100.715118”" 


9940-W 112.949112° 


FIGURE V-ic. Directions of Maximm Increase (Gradient) of TD's at R-7. 
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9940-W 7.762177° 


9940-K 344.397217° 
9940-M 38.852295° 


9940-Y 101.604533° 


FIGURE V-2a. Azimuths From Luces Point to 9940-M,W.X,ard Y. 


9940-Y 250.228424° 


9940-W 203.307236° / 


f 9940-X 191.624756° 


FIGURE V-2b. Directions of Hyperbolas at Luces Point. 


9940-Y 340.228424° 


9940-X 101.624756° 


9940-W 113.307236° 


FIGURE V-2c. Directions of Maximm Increase (Gradient) of TD's at IP. 
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2. Lanewidths 


Lanewidths were calculated from Equation IV-3a using the 
aZimuths computed in Section V-A-1. Note that the units of lanewidth 
are m/usec. With a standard lane of 1 usec the width of a lane is the 
same number in meters. Velocity of transmission was calculated using c 


and the index of refraction specified in the table of Loran-C constants 


supplied by the Defense Mapping Agency (Ref. 4). Sample calculations of 


Tanewidth follow: 


C 
— (V-1] 
n 


Velocity of Transmission 


2.99792453 x 108 m/sec 
1.000338 


= 2.996911624 x 108 m/sec 
Lanewidth of 9940-Y at Range-7 for standard 1 usec lane: 
(107~©sec) (2.996911624 x 108 m/sec) 
101.793397° - 38.358551° 


2sin 


285.024 m 


Results of the calculation of lanewidths are given as the L 
values in Table V-3. In the table, Iy, ly, amd Ly, correspond to the 


lanewidths for 9940-W, X, and Y, respectively. 
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3. Average TD's, Standard Deviations, Variances 


Because of the large amount of data involved, it was necessary 
to compute statistical information using a computer program. Since 
data was in the form of three TD's (one for each master-secondary pair) 
every five seconds, amd readings were taken for periods of 20 to 40 
minutes at each station, a Fortran 77 program with subroutines for 
statistical information was written (see Appendix C). 

Station Range-7 was the "known" station of the differential 
system, at which more readings were taken than any other. For this 
reason, and because regression equations of TD vs Time were needed only 
at that station, statistical information was calculated only for 
station Range-7. 


Mean of the TD's was computed by the following (Ref. 18): 


ZTD} 
Mean TD = [V-2 ] 
n 





Variances of the TD's were calculated using the following (Ref. 18): 


STD,“ - (ZTD;)2/N 
soo [V-3] 
N-1 


Standard deviations were calculated from variances by (Ref. 18): 


formula (Ref.18); 


Ss = (s“)% [v-4] 
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Equations V-2, V-3, and V-4 were incorporated into a Fortran 77 
computer program which was used to calculate statistics of TD's for 
9940-W, X, and Y. 

It should be noted that these statistics were in terms of time 
differences, while in some cases it was necessary to know standard 
deviations and variances in terms of distances. Equations IV-5 were 
used for this conversion where necessary. 

Corrected vs Uncorrected Data: One use of the statistics of 
the data was to compare the standard deviation of raw data to the 
standard deviation of data which was corrected for three secondary 
stations. It was considered unlikely that the three station correction 
could remove bias from the TD readings (error in terms of absolute 
position) because most of this is due to systematic error rather than 
random variations. It was believed that the three station correction ° 
would be most effective in reducing the effect of random variations in 
time differences. The statistics on uncorrected and corrected data show 


results consistent with this expectation (Table V-7). 


4. Error Ellipse am Circular Error at Range-7 


Range-7 was used as the known point of the differential system, 
and a larger amount of data was collected at that point. Using the 
crossing angles, standard deviations and variances of 9940-X and Y (the 
best pair at that point), dimensions and orientation of the standard 
error ellipse and CEP 50% and 90% circle radii were calculated for the 


point using the information outlined in section III-B. 
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The standard error ellipse and error circles for 50% and 90% 
were plotted over a scatter plot of data from Range-7. The result is 
shown in Figure V-3. 

Calculations of the standard error ellipse dimensions were as 


follows: 


Using 9940-X,Y: 


01 = Ogg49-y = (0.072546) (285.024) = 20.677 m 
O> = Ogg4p-x = (0.068638) (322.024) = 22.103 m 
@ = 250.075974° - 190.715118° = 59.360856° 


Calculating the semi-major and semi-minor axes of the ellipse: 


20.6772+22.1032+( (20.6772+22. 1032) 2-4 (20.677222.1032sin2¢) ) 7 
. pennies eee ae ee 
x 


2sin-¢ 
Oy = 30.595 m 
: 20.6772+22.1032-( (20.6772+22. 1032) 2-4 (20.677222.1032sin2¢) ) % 
0. eee — ee 
Y 
2sin-o 
Oy = 17.362 m 


Orientation of the major axis of the ellipse: 


20.677sin 2(59.360856°) 
tan 26 = 
20.677" cos 2(59.360856°) + 22.1032 


8 = 26.474268° 


Azimuth of major axis = 250.075974° - 26.474268 = 223.601716° 


223.601716° - 180° = 43.601716° (from North) 
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Figure V-3. Scatter Diagram, Standard Error Ellipse, and CEPsg amd 
CEP99 Error Circles, for Range-7 Data (see text for details). 


a7 


Using tables in Reference 6, pp.18-20, radii of the 50% ard 


90% error circles were calculated as: 


CEPS 9 = 27.936 m CEPo99 = 63.226 m 


Because the data included in the plot of Figure V-3 includes 
two intervals of data collection, the points are not as evenly 
distributed as they would otherwise be. For this reason there are 
fewer points near the center of the plot than would normally be the 
case. The error circles and ellipse seem to be reasonably correct for 


the given data. 


5. Correction Factor for Three Secondary Stations 


The correction for three secondary stations was designed as a 
technique to reduce the random error associated with time difference 
data. It was expected to increase the precision of the position rather 
than the absolute position. 

For a variety of reasons there will be some covariance between 
concurrent time difference readings from different master-secondary 
pairs. The correction for three secondary stations is not expected to 
reduce the error associated with this type of random error. The non- 
correlated random error in TD readings should be reduced by this 
correction to some extent. 

In order to reduce the number of computations and the time 
required to process data, most constant factors associated with the 
correction factor can be computed in advance and used as constants in 


the computer program. This reduces the correction to only a few lines 
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of relatively simple computation. Because of the large amount of data, 
this reduces significantly the time required to compute corrections. 

In the calculations of adjustments to positions of Range-7 and 
Iuces Point, the only factors which differ significantly are the long 
term time differences at the two points. In this case the sam Aj 
factors can be used for correcting both positions for three secondary 
stations. In cases where same of the parameters may differ 
Significantly, it would be necessary to calculate different A; factors 
for each position. Table V-3 gives parameters for input to the Fortran 
77 correction program for calculation of the A; parameters and 
calculation of the corrections to data. 


TABLE V-3. Parameters used in calculation of correction 
factorak. 


iiemo40112- 
100.715118° 
340.075974° 
0.150606 psec 
0.068638 psec 


0.072546 psec 


563.934 m/sec 


322.966 m/usec 


285.024 m/psec 
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6. Regression Equations for Time Difference vs Time 


Since only one Loran-C monitor receiver was available, it was 
necessary to develop a method of estimating TD's at the known point 
during times when differential corrections were required at the unknown 
points. Linear regression equations of TD vs time were used for this 
purpose. 

Using equations IV-21 and IV-22, the data from Range-7 was 
processed to obtain the required summations of time and time difference 
necessary to calculate the by and b, coefficients of the linear 
regression equation. Only Range-7 data from the nearest time periods 
before and after the differential correction period for Luces Point 
data were used. 

The procedure for processing the data was first to apply the 
correction for three secondary stations to the Range-7 and Luces Point 
data during the time periods of interest. Range-7 data was then 
processed to obtain the necessary summations for the calculation of the 
b parameters of the regression equations. Using the regression 
equations for Range-7 data and long term TD averages at the two points, 
differential corrections for the Luces Point data were calculated as 
described (Section V-7). Fortran subroutines (Appendix C) were added 
added to the statistics program to obtain the summations and a separate 
program was used for the differential correction of Luces Point data. 
Results of the corrections were compared to calculated values of TD's 
provided by the Defense Mapping Agency, Hydrographic/Topographic 
Center. TD vs time is plotted in Figures V-4 along with the regression 
line for the data. Regression equations are in Table V-4. 
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(usec) 


TD 


MiMi vss V8 
9940-W 9/28/87 


$6316.54 


$e316.4 





9.0 40.0 11.0 42.0 43.0 14.0 
TIME (PST) 


Figure V-4a. TD vs time and the regression line for 9940-W 
for Range-7 data corrected for three secondary stations. 
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Pye iic, Was ib 
9940-X 9/28/87 


GLI, 


Picce ae 


(ule ee) 


eS, 


Eysee.1 





27822.0 
8.0 10.0 41.0 42.0 13.0 14.0 


TIME (PST) 


Figure V-4b. TD vs time and the regression line for 9940-xX 
for Range-7 data corrected for three secondary stations. 
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(use 


eo 


TIME vs TD 
9940-Y 9/28/87 


42743 .6 


4c749.4 





42748.2 
9.0 410.0 411.0 42.0 13.0 44.0 


ete al) 


Figure V-4c. TD vs time and the regression line for 9940-yY 
for Range-7 data corrected for three secondary stations. 


Go 


TABLE V-4. Regression Equations of TD vs Time at Range-7. 


9940-W: TD, 16316.51564 - 0.00807243t 


9940-X: TD,, 2/522.38729 “=~ O.0IB S519 


9940-Y: TDy 42749.29371 + 0.01056707C 





time units: hours with minutes and seconds in decimal hours. 
7. Differential Corrections 


Equation IV-20 describes the differential correction which is 
applied to TD readings at an unknown point based on TD's at a known 
point. In this experiment the positions of all points are known, based 
on their published positions from survey data (see appendices). One 
point, Range-7, was chosen as the known point for purposes of 
differential corrections. Data from one other point was chosen as the 
unknown point for purposes of comparison to its known position. 

Three things must be known to calculate a differential 
correction based on equation IV-20. Rather than use long term averages 
in this equation, the expected TD's at both points were computed by a 
standard Loran-C program. Current TD's at Range-7 were calculated from 
the regression equations described in the last section. The ASF 
corrections in equation IV-20 were calculated for each point by the 
Defense Mapping Agency using a computer program which is used in the 
preparation of ASF correction tables, but which can be used for 


determining the correction at individual points. The difference in ASF 
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corrections is used in equation IV-20. Table V-5 gives expected TD 


values (used as long term averages) at the two points of interest: 


TABLE V-5. Calculated TD's used as long term averages (in psec). 


Range-7 
16314.68 16300. 09 


9940-X | PA pe Xe 27496.79 
9940-Y 42749.49 42755.87 





ASF corrections were approximately the same at both points que 
to the nearness of the points. Differences in ASF corrections used in 
equation IV-20 were very small: 0.005, 0.002, and 0.003 psec for 9940- 
W, X, and Y, respectively. Corrections of this magnitude are probably 
less than the accuracy of the data used to calculate them, but at 
longer ranges these may be large enough to be significant. 

The final equation for the differential corrections based on 
TD's at Range-7 are in Table V-6. Note that t is Pacific Standard Time 


in hours, with minutes and seconds expressed as a decimal. 


TABLE V-6. Differential Corrections (Cg in psec, t in hours PST) 


-1.83064 + 0.00807243t 


= ©0:79575 + ©6.01385199t 


OE19329 = 0.01056 707e 
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8. Accuracy of Position Fix 


Figure V-5a shows a plot of uncorrected Range-7 data based on 
expected TD's calculated from the known geographic position. The grid 
is in meters, oriented north-south and east-west, and with the center 
at the published geographic position of Range-7. When the Correction 
For Three Secondary Stations is applied and the average of the Range-7 
data is placed at the origin, the result is as shown in Figure V-5b. 

Note that the plot is now centered on the origin due to the 
fact that all absolute error (bias) has been removed from the data. 
The three station correction has also decreased the standard deviation 
and variance of the Range~-7 data as shown in Table V-7. 

Figure V-6a shows LIuces Point data before application of any 
corrections. These points are plotted based on the expected TD's at 
the known coordinates of the point. The center of the uncorrected 
Luces Point data is at x = 375.820 m, y = 85.501 m, and at a range of 
385.423 m and at an azimuth of 77.183° (from North) from the origin. 
Radii of the CEPsg and CEPgg error circles are 14.9 m and 28.4 n, 
respectively. 

Figure V-é6b shows Luces Point data after the application of 
both three station and differential corrections. Again, the origin is 
at the known position of the point. For visual comparison purposes, 
the same scale is used for both Figures V-6a and b. For the corrected 
data, the center of the data is at x = 47.213 mand y = -7.015 m, with 
a range of 47.731 m and at an azimuth of 98.451° (from North) from the 


Origin: CEPs5g and CEP9g radii were 12.6 m and 26.0 m, respectively. 
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UNCORRECTED 9940-x, Y 
— SCATTER DIAGRAM 
STATION RANGE-7 


200 





SCALE IN METERS 


Figure V-5a. Scatter Diagram For Range-7 Uncorrected Data. Scale is 
in meters oriented N-S and E-W, centered on the time 


differences calculated from the published coordinates 
of the point. 
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CORRECTED 9940-~x, Y 
300 SCATTER DIAGRAM 
STATION RANGE-7 


200 





ee 100 200 300 
f Ce SCALE IN METERS 


Figure V-5b. Scatter Diagram for Range-7 Data Corrected for Three 
Secondary Stations and Differential. Scale is in 
meters oriented N-S and E-W. Center of data is 
corrected to correspond with calculated time 
differences for the published coordinates of the point. 
See text for details. 


68 


UNCORRECTED 9940-xX, Y 
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Figure V-6a. Scatter Diagram For Luces Point Uncorrected Data. Scale 
is meters oriented N-S and E-W, centered on the time 
differences calculated from the published coordinates 
of the point. 
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CORR = Gil] bo S4 Oa 
SCATTER DIAGRAM 
STATION CBGes eat 


SCALE IN METERS 





Figure V-—6b. Scatter Diagram for Luces Point Data With Differential 
Corrections Based on Range-7 Data and With Correction 
For Three Secondary Stations. Origin of plot is at the 
time differences calaqulated fram the published 
coordinates of the position. 
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TABLE V-7. Standard Deviation and Variance of Range-7 Data With and 
Without the Correction For Three Secondary Stations. 


ie | tee oem | oertin 


Oo 0.068638 0.051161 





9940-X 
o2 0.004711 0.002617 







0.072546 0.052081 






0.002712 





0.005263 





B. Conclusions 


Based on the results given in part A, the absolute accuracy of 
the Luces Point Data was improved from an average error of 385.423 
meters to an average error of 47.731 meters. Practically all of this 
improvement was a result of application of the differential correction. 

Precision of the Luces Point data was improved from 14.9 m and 
28.4 m for the 50% and 90% CEP radii, uncorrected, to 12.6 m and 26.0 m 
for the corresponding radii, corrected. Most of this improvement is a 
result of the correction for three secondary stations. 

These figures are based on the processing of 486 sets of TD 
data from Range-7, the "known" point of the differential system, and 
252 sets of TD data at Iuces Point. Each set of data included time and 
three TD's, one for each master-secondary pair (9940-W, X, amd Y). 

Under ideal conditions a better absolute correction would be 
expected for these points. From an examination of the characteristics 


of the data, the weakest point is probably the lack of a continuous set 


fan 


of TD readings at Range-7. Error in the estimate of TD's based on the 
linear regression lines may be as much, or more than, the 48-meter 
average error in the corrected positions calculated for Luces Point. 

Other sources of error may include the fact that some of the 
Loran-C signals traversed some seawater path near Luces Point, with 
unknown results on TD readings at that point. This type of error is 
known to result under similar conditions, but no means of estimating 
its magnitude is available at this time. 

In summary, the absolute positions calculated from differential 
corrections of Loran-C data seem to be a significant improvement over 
the positions calculated from raw data. While the aaoeseer for three 
secondary stations does not add much to the estimate of positions in 
this experiment, the improvement of 50% CEP radii by about 15% does 


seem Significant enough to warrant further consideration. 
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APPENDIX (C-|] 


THIS PROGRAM READS DATA FROM A FILE OF TIME, A, B, C 
CONTINUES TO THE EWLD OF THE FILE, CALCULATING THE MEAN, 
VARIANCE, AND STANDARD OCEVIATION OF THE A, B, AND C COLURNS, 
AND PRINTS OUT THE RESULTS INTO A FILE todout.f 
Gouble precision tiwe,timel, tiwe2,tine3,a,b,¢,nm,na,nb,nc, 
*SuB,$uUBG, SUBD, $UBC, SU@a2, SUBb?2, SUBC?2, SUBAl, SUBD}, Sueci],al,bi,c) 
open (unitei,file = ‘out9.t', form = ‘formatted’, 
eaccess: sequential’) 
open (units2,file = ‘out6.¢’, fore = ‘formatted’, 
saccess= sequential’) 
revind 1 
react) .°) {rstin 
reed(1,*) secIn 
write (*,4) timei,time2,time3,a01,61,¢1,nm,nma,nb,ne 
time = timel + time2/60. «¢ time3/3600. 
write (2,8) time,el,bl,¢l 
Sum © sue ¢ }. 


sume = sume + 6 
sume) = supal ¢ a] 
suma2 = suma2? + ae 
su@eb = sumb + b 
Sumb) = sumti + bil 
Sueb?2 = sumbd2 ¢ b"*b 
SuRC = gumc ¢ € 
sume) = sume] ¢ c) 
susc2 © susc2 ¢ cc 


read(i1,2,errz=lO0,iostatzioerr) timel,time2?,time?,ei1,a,b61,D, 
+i Colne 

it (yoerr? 16,50, 70 

Sur = $um - j. 

write(*®,6) joerr 
call s2(suma2?.sumea,sun,sumai) 
coll s2{sumb2, sub, sum, sumb} ) 
call s2{sumc2.suec, sum, sueci} 

formatitioe?. O..téend..0,t7,02-0, the, de. ceteteGe vote 7, 
d8.3,t030,05.3,056.0855.tS9, 05.5) 

formatt (5.0.27. -0, 2s, 45,0, 2), 1 Oss ca lel yn Ons. ge 
416,46,16,i16) 

format(’fioerr « °,i3) 
Yorsestit)2.6,.2%,110.3,2:,110.5,28. 1 tos) 

format(“date error’) 

stoo 

end 


+ 


Subroutine s2(12,1,4,$81) 
double precision 12,1,6,@ean,vaerian,sd,$] 
meen = si/a 
varian = (12-2%1/0)/{6-1) 
so = dsart{varian) 
write(*,12) wean, varian,sad 
write(2,12) mean, varian,sa 
return 
format(’sean = °,f12.3,21, variance = °,f12.8,21,'s0 = ',712.8) 
end 
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APPENDIX C-2 


ree ---ce------- oe ewe ween eee eee eee See ee ee eee ee eee 
rea --- PROGRAM FOR GAUSSIAN ELIMINATION WITH PIVOTAL CONDENSATION ------ 
FeB -- errr rete FOR SOLUTION OF SYSTEMS OF LINEAR EQUATIONS -------------- 
ae WiTH UP 710 10 UNKNOWNS = = one cece mre enn nen ere cern 
7) oe ee ee 
res 

open “i°, @1, “ytest.f° : dim AQ(10,11):die Baelii):die Me(iididie xreldl) 


input @1, rows!, cols! 

ros fous! = ¢ s cols! =: roul = 1 4 e€oll = 1: highrow «= 1 
ren 

rem <-er-- input the eatrixz: 

ree 

for j = 1 to rows! 

for k «© 1 to cols! 

input @1, A&(j,k) 

if k () cols! then 160 

Ag(j,k) © -1 * aalj,k) 


Nnert k 

next j 

ree 

Gem —--=== BEGIN REIN PrOOr Bw === <= <9 < ea 6a eRe a a= Ae Sam = = = nim = 
rea 

gosub $00 : rem --- check if uniaue solution 

gosub 890 : rem --- prints out erreay 

gosub 680 : reg --- checks for row with highest first non-zero entry 
ree 

if highrow = rowl then 260 

rewr 

gosut 700 : ren --- erchenge rows if necessary 

gosut 760 : rem --- sets first nom zero entry in rowl to 1.0 
gosub 820 : rem --- sets first entry to zero in tollowing rows 
ren 

rowl = rowl ¢ 1: ¢oldl s coli1 + 1 : highrow = rowl : rem --- Increment counters 
ree 

if rowl () rows then 200 : rem --- checks if lest row hes been reached 
gosub 760 : rep --- sets first non zero entry in final roe to 1.0 
gosub 890 : rep --- fine} printout of array 

rem 

rem e+e-- calculates the unknown 1: velues: 

ree 


r@(cols!) «= 1: o = rows! 
for : =» @ to il step -1 


for 3 = 4 to cols! - 1 

yeCi) © sali) - sa(je3) ° aali, jet) 

nert j 

nert {i 

ree 

rem e+e prints out the 3: values: 

rem 

tor k = 1 to rows! 

print “B° kel ° « ° xefk) 

next k 

end 

rep ----- end of main program ---- mmm mr mr rrr rr rr rr rr eee ene 
res 

rem cece begin subroutines onc - cen nce eee nee eer nee renee 
ree 


Tor 4 = 1 to rows! 

for 3 = 1 to rous! 

if 1 = j then 610 

Tor eos ) to cols! 

ff as(i,k) « O then 560 
Me(k)= as(j,k)/as(1,k) 
next k 
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for wo © 1 to cois! - 1 

if m@(a) () Ae(mel) then 610 

neat @ 

print “EQUATIONS NOT JNOEPENDENT ... WO UNIQUE SOLUTION ... OWE SOLUTION FOLLOWS’ 
next j 

neit i 

return 

ree 

rem ---r- subroutine to check for row with highest first non-zero entry: 
for i = rowl to rows! - {i 

if abs(ael(i,coll)) ») ebs(atliel,coll)) then 680 

highrow = ie 4} 

next i 

return 

rem ----- sunroutine to exchange rows and put highest first entry first: 
for { = coll to cols! 

Beli) = A@frowvi, i) 

Atlrowi,i) = at(nighrow,i) 

At(hionhrow, i) = Beli) 

next i 

return 

re@ ----- subroutine to set first entry in first row to 1.0: 
R@ = a@(rowi,coli) 

for k = coii to cols! 

if a@(rowi,coli) = 0 then 310 

AS(rowi.k) = a@(rowl, k)/R8 

print bk 

neat k 

return 

rem <----- subroutine to set other first entries to zerec: 

for j = rowl « 1 to rows! 

Fe = AB(j, coll) /AB(rovl,coll) 


for k = coll to cols! 

AS{j,b) = AB( j,k) - ARLroul,k) PS 

nert k 

neat j 

return 

rem ----+- subroutine to orint out the matrix A: 


for Jj © 1 to rows! 

for k © 1 to cols! 

print ae(Jj,k), 

next k 

print 

nert j 

Print POW); “enn eee mee ee nee ene ee ne eee ee ee en rn nn enn rn re ren re ren eee- 
return 

print “NOT ENOUGH EQUATIONS ... MANY SOLUTIONS ARE POSSIBLE” 
return 

end 
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PROGRAM aAcr 
A= €3782C66 6 
B = 6*Sése:.8 
OFEN(UNIT=1. FILES azin. f°, FORM: FORMATTED’, 
*aACCESS* SEQUENTIAL’) 
REWIND j 
Beimutlr is from point 1 to point 2 
READ(1,2) ALATZ, ALON? 
READ(1,2) ALAT1,ALON1 
feaGc = 0.0 
CALL JNVEA,B,ALAT1, ALON] ,ALAT2,ALONZ, FLAT, DIST. AZ) 
AZ = Ag « 3.7159 2¢ Sz 
AZ = AZ/2.14358926$6°1EC.0 
IF {A2-360.0} 16,310.20 
AZ = AZ - 360.0 
WRITE SoS) AZ 
m= ALAT 2) © 169592¢56"° 760. 0 
Y ALONT/3.141592¢86°189.0 
Voz ALATISS. yoUS ere foyer 
Ue akon Solel SS2e 50" 126) f 
WRITE!*, 6) XY 
BETTER". S) Vu 


ZOE AUES Sela 9, £1 


163 


iC: 


VO’ 
Dee 


Beaker) CPSTANCE = 2f 2 Noo METERS. Sa SAS SALT: © oP ee oe 


7 ENTE TR) ) 


Riemer ort TNS sek EAT OO wo Pore.§. No Sek SOLON, 20 Peo 
POGmene tr CLs the. wee LATS Dy eS ose LGN: | a: Qi kG 
EN co 


SUGEGCHI ThE IhvlE OFAC FORRG BLT. FLNYOREYISFEN, <CFRLG.EIStaee] 
Pete oet To4 hei. 1s 

FLA? © 1 G-FORAC EGF AL ° 
FEELS = PLETED ET 

ae FAUATC (T° 25 

7 Ap ae eae Wee area 

Pe Tay See ae 
ae oo lala Smo ala 

Pos FAAS "O. O¢: © 

Fe SIA ee ery 

Fo  latacs| Fa 

Ripper hetertg aoc 

Pie 3. 1678S 2€ $2 

VW ieee Gage Sees: 

CEL Gs. 1". 6 

PEVAIe= ATARIS 1. O-FUSTICSINIF C TINSCOSURL Ta) ) 
SEETA} = SINICETA)) 

CEETA] = COS(EETA]) 

FETA? AARC IGF UAT) “SINIFLT2)/COSIRET?)) 
SEETAD = SIN(BETAS) 

CEETAZ = COSI BETAZ) 

CELL = FLNI-FLEG 

ADEUL = AESITELU) 

VESRPAGE bie F iolmeroe TCE. IGE 

ACEUL = THWOFT=ACELL 

SIUEL SIN ASELL) 

CODEL = COSTADELL 

(Yee itis Sah Silas ie 

CREA! CEST Kz 

HIT = A+E*CODEL 

H 


€ 


& PRGSUbE ChE Th.) (29 (SEE TA2 CEE TAI=SEETAI“CEETAs COUEL I> * 2 
See S)GrLy SIP HI 

Peas > C 

HJ] = ASITNISJPH]) 
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213 
219 


IF (COPHI) 207,208,208 
PH] « PI-PH] 

PHISO = PHI *PHI 

CSPH] = 1.07513FH1 
CYPHI © COPHI/SIPHI 
PSYCO = SIPHI]*COFHI 


APPENDIX (C(-3 


EN"ER"(PS°CPHIVPESYCO)+F2*FHISO*CTIPHI-Fé"PSYCOUCCDEHL cocna 


TERM) = F7°PHI 

TERM2 & AT(FE*SIPHI-F2*PHISO*CSPH] ) 

TERMS & EM*(F2*°PHISO*CTPHI-FB* (PHI +PSYCO) ) 
TERMS = ATA*F2°*PSYCO 

TERMS «& 

TERME = ATEM*F2*(PHISOSCSPHI+PSYCO*COPH! ) 


DIST £& PORAD*(TERMI*4TERM2+¢TERMI-TERMACTERMS*TERMO ) 


TERM? & FE*PHI 


TERMB © A*(F2*SIPHIOFALAT2°PHISO*CSPH] ) 


TERMD & EM*(F3*FS*PSYCOCFALAT2°PHISO*CIPHI-FL PHI ) 


ZLAM = C*(JERM7-TERMB*TERMS) ADELL 


CTAZ=(SBETA2*CBETAI-COS(ZLAMZ)*SBETAI*CBETA2)/(SIN(ZLAM)*CBETA2) 


IF (CY¥AZ) 210,209,210 
ctaZ = 0.00000005 

AZ = ATAN(1.0/CTAZ) 

IF (OCU e2ISy2be, 216 

if (DEUL-F1) 21i2 212, 212 
TF UCTAZ) 22052214224 

AZ = AZ+*P] 

GOTO 221 

LE UECYAZ) 217,278. 216 

FE (OECD Lee eee 
TF CCTAZY 212,208 -218 

AZ = PI-A2 

GOTO 271 

AZ = TWOPI-A2 

AZ = AZ<¢F] 

AZ = AZ-TWOP] 

BE UA 23219219 

AZ = AZeTWOFP] 

RETURN 

END 
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